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ABSTRACT 

Recent observations have revealed that Type la supernovae (SNe la) are not perfect standard candles 
but show some variations in their absolute magnitudes, light curve shapes, and spectra. The C/0 
ratio in the SNe la progenitors (C-0 white dwarfs) may be related to this variation. In this work, we 
systematically investigate the effects of stellar mass (M) and metallicity {Z) on the C/0 ratio and its 
distribution in the C-0 white dwarfs by calculating stellar evolution from the main-sequence through 
the end of the second dredge-up for M = 3 — 9Mq and Z = 0.001 — 0.03. We find that the total carbon 
mass fraction just before SN la explosion varies in the range 0.36 ~ 0.5. We also calculate the metallicity 
dependence of the main-sequence-mass range of the SN la progenitor white dwarfs. Our results show 
that the maximum main-sequence mass to form C-0 white dwarfs decreases significantly toward lower 
metallicity, and the number of SN la progenitors may be underestimated if metallicity effect is neglected. 
We discuss the implications of these results on the variation of SNe la, determination of cosmological 
parameters, luminosity function of white dwarfs, and the galactic chemical evolution. 

Subject headings: stars: evolution — stars: interiors — stars: white dwarfs — stars: supernovae: 
general — stars: abundances — nuclear reactions 



1. INTRODUCTION 

Type la supernovae (SNe la) are characterized by the 
absence of hydrogen and the presence of strong Si II line 

■ in the early time spectra. Majority of SNe la have almost 

■ the same absolute brightness, thereby being good distance 
, indicators for determining cosmological parameters (e.g., 

Branch and Tammann 1992; Sandage and Tammann 1993; 
Nomoto, Iwamoto, & Kishimoto 1997 for a recent review). 
However, SNe la are not perfect standard candles but show 
some variations in their absolute magnitudes, light curve 
shapes, and spectra (Phillips 1993; Branch, Romanishin, 
& Baron 1996; Hamuy et al. 1996). It is important to 
understand the source of these variations and the depen- 
dence on the metallicity in order to use high-redshift SNe 
la for determining the cosmological parameters. 

It is widely accepted that SNe la are thermonuclear ex- 
plosions of accreting C-0 white dwarfs (WDs), although 
detailed explosion mechanism is still under debate (e.g. 
Ruiz-Lapuente, Canal, & Isern 1997). Previous studies 
have shown that typical SNe la are well described by the 
Chandrasekhar mass model, in which accreting C-0 WDs 
explode when the masses of these WDs approach the crit- 
ical mass Mia — 1-37 — 1.38M0 (Nomoto, Thielemann, 
& Yokoi 1984; Hoflich & Khokhlov 1996). The explo- 
sion is induced by the carbon ignition in the central region 
and the subsequent propagation of the burning front. Be- 



hind the burning front, explosive burning produces a large 
amount of ^^Ni. There are two models for the burning 
front propagation: the deflagration model (e.g., Ivanova 
et al. 1974; Nomoto, Sugimoto, & Neo 1976; Nomoto et 
al. 1984) and the delayed (or late) detonation (DD) model 
(Khokhlov 1991; Yamaoka et al. 1992; Woosley & Weaver 
1994). 

There are several possible sources of the variations of 
brightness and spectra: for example, variations of the 
age of the WDs and the mass accretion rates, which af- 
fect the ignition density, variations of the WD progen- 
itor's main-sequence mass and metallicity, which affect 
the C/0 ratio in the WDs, and the variation of the 
deflagration-detonation transition density (ptr) in the DD 
model. These effects modify the explosion energy and the 
synthesized ^^Ni mass, which affects the maximum bright- 
ness and shape of the SN la light curve. 

Regarding the C/0 ratio, previous explosion models 
have been constructed by assuming a uniform mass frac- 
tion ratio C/0 = 1, though previous results indicated C/0 
< 1 in the central region (e.g. Mazzitelli, & DAntona 
1986 a,b; Schaller et al. 1992, SSMM hereafter). This is 
partly for simplicity and partly due to its strong sensitiv- 
ity on the treatment of convection and the ^^C(a,7)^^0 
rate. However, the effect of the different C/0 ratio is 
not negligible. For example, recently Hoflich, Wheeler, & 
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2 Evolution of 3-9 M© Stars from Z=0.001 to 0.03 

Thielemann (1998) constructed explosion models by arbi- 
trarily changing the C/0 ratio to C/0 = 2/3 and found 
some significant effects on the light curve and spectra. 

Intermediate mass stars by definition burn only hydro- 
gen and helium, and form C-0 cores after He burning. 
These cores will become C-0 WDs, thus becoming SNe 
la progenitors in certain binary systems. Previously, sev- 
eral groups have investigated intermediate mass stars, but 
they have not paid much attention to the C/0 ratio and 
the C-0 core mass size with its systematic dependence on 
the stellar mass and metallicity, as well as its effects on 
SNe la. For example, SSMM and a series of subsequent 
papers (Schaerer et al. 1993; Charbonnel et al. 1993) 
studied the influence of the metallicity on the life times, 
luminosities and the chemical structure. However, they 
only showed the central and surface chemical abundances 
and did not show the C-0 core mass sizes. Salaris et al. 
(1997) showed the radial profile of the C/0 ratio in the 
intermediate mass stars as well as the C-0 core mass size, 
but only for solar metallicity. Vassiliadis & Wood (1993) 
calculated C-0 core mass sizes for the stars below 5 Mq, 
but did not pay much attention to the C/0 ratio. The 
works prior to 1992, such as Castellani et al. (1990), Tor- 
nambe & ChiefR (1986), Becker & Iben (1979), used Los 
Alamos (Huebner et al. 1977) or Cox & Stewart (1970) 
opacity tables, and hence their results would be signifi- 
cantly different if latest OPAL opacity (Iglesias & Rogers 
1993) was used. 

In this paper, using updated input physics, we evolve 
the intermediate mass (3 — 9Mq) stars and examine the 
stellar metallicity and mass dependence of the C/0 ra- 
tio, its radial profile, and the C-0 core mass size. These 
quantities arc sensitive to the treatment of convcctive mix- 
ing; however, for the same treatment of convection, the 
^^C(q!,7)-^^0 rate dominantly determines these. As de- 
scribed in section 2, we calibrate this rate by the massive 
star nucleosynthesis. 

We also obtain metallicity dependence of zero-age main- 
sequence (ZAMS) mass ranges of SN la progenitors, and 
especially its upper limit Myp, where Myp is the ZAMS 
mass at the border between the star that ignites carbon 
under non-degenerate condition and the star that forms a 
strong degenerate C-0 core. This quantity is important 
for various applications, such as galactic chemical evolu- 
tion and the luminosity function of WDs in our galaxy. 
However, earlier determination of A^up was based on old 
radiative opacity tables (e.g. Becker & Iben 1979, BI79 
hereafter; Castellani et al. 1985, CCPT hereafter; Tor- 
nambe & Chieffi 1986; Castellani et al. 1990) or not by 
exact calculations (Schaerer et al. 1993). Also the cal- 
culations during the 1980s and early 1990s often assumed 
strong overshooting or breathing pulses, which may not be 
preferable with the use of the OPAL opacity. Their results 
are qualitatively consistent with our results but quantita- 
tively different. Our Mup is significantly smaller than the 
estimates in BI79 for all Z we calculated, and in Schaerer 
et al. (1993) for Z< 0.01, and it shows a large decrease 
toward lower metallicity. Such a dependence would be 
important for the distributions of SNe la brightness as a 
function of metallicity. 

This paper is organized as follows. In the next section, 
our numerical code, input physics and parameters are de- 
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scribed. Section 3 shows our results and compares with 
previous works. Their implications are discussed in Sec- 
tion 4. 

2. NUMERICAL CODE AND ASSUMPTIONS 

We modify and update the Henyey type stellar evolution 
code (Saio & Nomoto 1988). Important improvements in- 
clude the use of the OPAL radiative opacity (Iglesias & 
Rogers 1993), neutrino cmissivity (Itoh, Nishikawa, & Ko- 
hyama 1996) and a nuclear reaction network with larger 
number of isotopes. Most of the nuclear reaction rates are 
provided by Thielemann (1995). This network is coupled 
to the dynamical equations for stellar evolution implic- 
itly in a similar manner as SSMM and Weaver & Wooslcy 
(1993). That is, the network with a reduced number of 
isotopes is run during Henyey relaxation. Nuclear energy 
generation rates between the reduced and full networks 
are almost identical. After dynamics is converged, nucle- 
osynthesis is calculated with a full network code for He 
burning. The isotopes included are shown in Table 1. For 
H burning, the network code of Saio & Nomoto (1988) is 
used with updated nuclear reaction rates. 

Treatment of convection is one of the serious uncer- 
tainties in the theory of stellar evolution. We adopt the 
Schwarzschild criterion for convective instability and the 
diffusive treatment for the convective mixing (Spruit 1992; 
Saio & Nomoto 1998). This formalism contains a parame- 
ter /fc of order unity which is proportional to the diffusion 
coefficient, and also to the square root of the ratio of "so- 
lute" to thermal diffusivity. We found that the results 
were not sensitive to this parameter at least in the range 
fk — 0.1 — 1. In this paper, we set fk = 0.3. We use 
the mixing length parameter a = 1.5 for the outermost 
super-adiabatic convection. Chieffi, Straniero, & Salaris 
(1995) found that a larger value (a ~2.25) fits observa- 
tions well with the use of OPAL opacity. However, as 
shown in the next section, in our code the evolutionary 
tracks in the H-R diagram are closer to those of SSMM 
for a = 1.5. Also as described in the next section, the 
core property, especially the central C/0 ratio, is insensi- 
tive to this parameter. Therefore, our conclusions in this 
paper are little affected by the choice of a different a. We 
neither artificially enhance nor suppress convective regions 
from the Schwarzschild criterion throughout the calcula- 
tions, and hence we do not include convective overshooting 
and breathing pulses at the end of He burning. It is con- 
sidered that with the use of the OPAL opacity these effects 
are not so large as considered before (e.g. Stothers & Chin 
1991); however, whether these effects are significant or not 
is still uncertain. 

Once the treatment of convection is fixed, the 
^'^C{a, 7)^®0 rate dominantly determines the C/0 ratio af- 
ter He burning and the C-0 core mass size. Since the track 
in the H-R diagram is not sensitive to the ^^C(q:,7)^^0 
rate, it is very difficult to constrain this rate from the sur- 
face property of the stars, such as spectral observations of 
clusters. In this paper this rate is chosen to be 1.5 times 
the value of Caughlan & Fowler (1988, CF88 hereafter). 
This choice, denoted by 1.5x CF88, with the same treat- 
ment for the convection yields the central carbon mass 
fraction X{^'^C) ~ 0.20 after the He burning of the 25 
M0 star (Umeda & Nomoto 1998). Therefore, it approxi- 
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matcly reproduces the central C/0 ratio of the Nomoto & 
Hashimoto (1988, NH88 hereafter) model, being expected 
to be consistent with solar abundances of Ne/0 and Mg/ O 
in massive star yield (Hashimoto, Nomoto, & Shigeyama 
1989; Thielemann, Nomoto, & Hashimoto 1996; Wooslcy 
& Weaver 1995). Wc note that although the ^'^C{a,j)^^d 
rate is constrained most stringently by the nucleosynthesis 
argument of massive star yield (Weaver & Woosley 1993; 
Thielemann et al. 1996), our knowledge of this rate is still 
limited, because there arc still uncertainties in the chem- 
ical evolution, massive star evolution, and the explosive 
nuclear burning models. In Section 3.3, we examine the 
dependence of this rate by applying the rates of 1.7, 1.5, 
and 1.0 times the rate of CF88 and the rate of Caughlan 
et al. (1985, CF85 hereafter). More detailed investiga- 
tion of the C/0 ratio in massive stars, and its dependence 
on ^^C(a, 7)^0 rate and treatments of convection will be 
given in Umeda & Nomoto (1998). 

We adopt the solar chemical composition by Anders & 
Grevcsse (1989), the solar He abundance 1^0 =0.27753 
and metallicity .^©=0.02 (Bahcall & Pinsonneault 1995). 
We assume that the initial (zero age main sequence or 
ZAMS) helium abundance depends on metallicity Z as 
Y{Z) = Yp + {AY/AZ)Z (for Z < 0.02), where AY/AZ 
is a constant and Yp is the primordial helium abundance. 
In this paper we adopt Yp = 0.247 (Tytler, Fan, & Buries 
1996), and thus AY/AZ = 1.5265. For Z > 0.02, we only 
show the case of (Z,y)=(0.03, 0.28). Different Y yields 
different stellar evolution. Especially, the star is brighter 
for a larger Y. However, as shown below the effects of 
changing Y is not very large. 

We also include an empirical mass loss rate (de Jager, 
Nieuwenhuijzen, & van der Hucht 1988) being scaled with 
metallicity as {Z/0.02)°-^ (Kudritzki et al. 1989). 

3. RESULTS 

We calculated the evolution of intermediate-mass (3 — 
9Mq) stars for metallicity Z=0.001 - 0.03. Figure 1 shows 
the evolution in the H-R diagram and the central density 
(pc) vs. central temperature (Tc) diagram of our 5 Mq 
(Z=0.02) model compared with SSMM. The luminosity of 
our model tends to be lower than SSMM and pc is higher 
for the same Tc. There are many differences in the in- 
puts, which may be the reasons for these differences: we 
adopt Y = 0.2775 for Z = 0.02, while SSMM assumed 
y = 0.3. Our choice of the i2c(a, 7)1^0 rate (1.5 x CF88) 
are smaller than CF85 ~ 2.3 x CF88 in SSMM. We use 
the larger nuclear reaction network. The treatment of con- 
vection is different. Especially they assumed overshooting 
from a convective core, which we do not. Figure 1 shows 
that the difference of Y alone does not explain the discrep- 
ancy. Also, we find that the tracks on the H-R diagrams 
are almost identical even though the ^^C(a, 7)^^0 rate is 
different. Therefore, probably the treatment of convec- 
tion is mainly responsible for the differences. Indeed, the 
tracks on the H-R diagrams of the no-overshooting models 
in Bressan et al. (1993) are similar to our results. 

In this figure we also show the effect of mixing length 
parameter a for the outermost super-adiabatic convection. 
This parameter changes the location in the H-R diagram 
after a star becomes a red giant. On the other hand, the 
central Pc - Tc track is almost independent of this param- 



eter. As a result, the central C/O ratio is almost identical 
between the a = 1.5 and 2.25 cases. Therefore, the choice 
of this parameter is not so important for the purpose of 
this paper. Chieffi et al. (1995) found that a = 2.25 fits 
observations well; however, the evolutionary tracks in the 
H-R diagram are closer to those of SSMM for a = 1.5 than 
2.25, and hence we adopt a = 1.5 in the followings. 

3.1. Metallicity effects on stellar evolution 

In the ranges of stellar masses and Z considered in this 
paper, the most important metallicity effect is that the 
radiative opacity is smaller for lower Z. Therefore, a star 
with lower Z is brighter, thus having a shorter lifetime 
than a star with the same mass but higher Z. In this sense, 
the effect of reducing metallicity for these stars is almost 
equivalent to increasing a stellar mass. This is verified in 
Figure 2, which shows the H-R diagrams of the 5Mq mod- 
els for different metallicities (upper panel) and the models 
with Z = 0.02 for different stellar masses (lower panel). 

For stars with larger masses and/or smaller Z, the lu- 
minosity is higher at the same evolutionary phase. With 
a higher nuclear energy generation rate, these stars have 
larger convective cores during H and He burning, thus 
forming larger He and C-0 cores. This property is de- 
picted in Figure 3(a) and (b). 

Figure 3(a) shows abundances in mass fraction in the in- 
ner core of the 7 Mq star for Z=0.001 at age f=48.45 Myr, 
with the central temperature and density log Tc (K)=8.58 
and log Pc (g cm~^) =6.82. In Figure 3(a), the He burning 
shell is located around the mass coordinate Mr = I.IMq, 
which indicates that the C-0 core of this model will ex- 
ceed 1.1 Mq. Since this is well above the critical mass for 
the off-center carbon ignition, this star will ignite carbon 
in the non-degenerate region and will not form a C-0 WD 
but form an 0-Ne-Mg core (Nomoto 1984). 

If the star is less massive or has higher metal abundance 
than the model in Figure 3(a), the C-0 core can be smaller 
and the carbon ignition is avoided as shown in Figure 3(b) 
for the 7 Mq model of Z=0.02 at age t=51.87 Myr (at 
the end of the second dredge- up) . The dredge-up is a phe- 
nomenon during which the surface convection is penetrat- 
ing inwards and the inner material is being brought up to 
the surface. The first dredge- up occurs before or around 
the He ignition and the second one occurs during the He 
shell burning stage for intermediate-mass stars. During 
the second dredge-up phase, the thick He shell which is 
seen in Figure 3(a) becomes thinner and thinner because 
of the C-0 core growth through He shell burning and the 
penetration of the surface convection. Therefore, as shown 
in Figure 3(b), the He shell becomes very thin (in mass co- 
ordinate) at the end of the second dredge-up phase. Chem- 
ical abundances at the end of the second dredge-up for the 
6 and 4:Mq stars (with Z = 0.02) are also shown in Figures 
4 and 5, respectively. 

3.2. ZAMS mass ranges of SN la progenitors 

The stars that do not ignite carbon will evolve into the 
asymptotic giant branch (AGB) phase. During this phase, 
if they are single stars, the C-0 core mass increases slightly 
due to H-He double shell burning until they lose almost 
all hydrogen-rich envelopes by wind mass loss. If they 
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are in close binary systems, their hydrogen-rich envelopes 
will be stripped off in much shorter time scale due to the 
Roche lobe overflow, and C-0 WDs are left. When the 
companion star evolves off the main-sequence, mass trans- 
fer occurs from the companion star to the WD. A certain 
class of these accreting WDs become SNe la depending on 
the accretion rate and the initial mass of the WD (e.g., 
Nomoto 1982; Nomoto & Kondo 1991) 

The WD mass at the onset of mass transfer from the 
companion star (denoted by Mwd,o) is critically impor- 
tant for the fate of the accreting WD. For Mwd,o as 
low as < 0.7 Mq, the WD mass is difficult to reach 
Mia ^ - 1.38M0. For higher Mwd,o, the WD can 
become SNe la in a certain binary system (Hachisu, Kato, 
& Nomoto 1996; Li & van den Heuvel 1996). However, 
the MwD,o depends on various uncertain factors, such as 
a binary separation, effects of recurring He shell flashes, 
and a mass loss rate during the AGB phase. In the follow- 
ing discussion, we assume that AfwD.o is not significantly 
larger than the C-0 core mass at the end of the second 
dredge-up phase, which we denote Mco- 

Figure 6 shows the metallicity dependence of the rela- 
tion between the ZAMS mass (Af,„s) and Mco- We note 
that stars with higher metallicity form smaller C-0 cores if 
they have the same Mms. Our Mms vs. AIqo relations are 
similar to Castellani et al. (1990) which used Los Alamos 
Opacity (Huebner et al. 1977) and assumed overshooting 
but suppressed breathing pulses at the end of He burning. 
They calculated Z=0.02, 0.006, 0.002 and Y=0.23, 0.27 
cases for 3, 4, 5, 7, 9 Mq stars. Our M^b vs. Mco rela- 
tions are also consistent with Vassiliadis & Wood (1993), 
which used OPAL opacity without overshooting. 

Figure 7 shows Afyp as a function of metallicity, which 
is compared with BI79 and CCPT. Our results in these fig- 
ures are quantitatively different from BI79, because BI79 
used old Cox & Stewart (1970) opacity. Our results are 
also different from CCPT's because they assumed rela- 
tively large overshooting. Although Castellani et al. 
(1990) did not estimate M\jp exactly, their results indi- 
cate that their Afup lies closer to ours than those of BI79 
and CCPT. CCPT and Castellani et al. (1990) assumed 
strong overshooting, but we note that when OPAL opac- 
ity is used the overshooting effect should not be very large 
to meet the observational constraints (Stothers & Chin 
1991). Schaerer et al. (1993) also estimated Afup by com- 
paring their grids of stellar models with Tc vs. pc tracks 
in Maeder & Meynet (1989), which is shown in Figure 7 
by the crosses. Our results show a stronger metallicity de- 
pendence than Schaerer et al. (1993). Especially, their 
Afup at Z=0.001 is significantly larger than ours. 

Tornambe & Chieffi (1986) studied Afup for very ex- 
tremely metal-deficient stars (Z < 10~^) and found that 
the A^up increases for smaller Z. Though verifying this re- 
sult may be interesting, we do not treat these extremely 
low metal stars because stars with Z < 10"'^ may not be- 
come SNe la progenitors in the Chandrasekhar mass model 
(Kobayashi et al. 1998). 

In Figure 7, we also plot A/lOj which is defined to be the 
lower ZAMS mass limit for stars to become SNe la. Here 
we assume that in order for the accreting WD to reach 
Afia, Mco should be larger than a certain critical mass, 
typically ^ O.8Af0 (Hachisu et al. 1996). Since this crit- 
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ical Mco depends on the binary parameters, three cases 
of Mco are shown in Figure 7. 

In Figure H we show the number of SNe la progenitors, 
i.e. the number of stars in the range ATlo ^ M < Afyp, 
normalized to the case oi Z — 0.02 and Mco — 0.8Mq. 
Since A/lo is sensitive to the binary parameters, we con- 
sider four cases of A/lo in Figure 8. In obtaining this 
figure, Salpter IMF is used. This figure shows that if 
Salpeter IMF holds also in a lower metallicity environ- 
ment, the number of SNe la progenitors increase toward 
low metallicity, because stars with lower metallicity have 
larger C-0 cores for the same Af„is- It should be noted 
that these numbers do not necessarily proportional to the 
SNe la rate, because we need to take into account the bi- 
nary companion's lifetime in order to calculate the SNe la 
rate. 

3.3. C/0 ratio in white dwarfs 

From Figures 3-5, we find that the central part of these 
stars is oxygen- rich. The C/0 ratio is nearly constant in 
the innermost region, which was a convective core during 
He burning. Outside this homogeneous region, where the 
C-0 layer grows due to He shell burning, the C/0 ratio 
increases up to C/0 > 1; thus the oxygen- rich core is sur- 
rounded by a shell with C/0 > 1. In fact this is a generic 
feature in all models we calculated. The C/0 ratio in the 
shell is C/0 ~ 1 for the star as massive as ~ 7AfQ (Fig. 
3b), and C/0 > 1 for less massive stars (Figs. 4 and 5). 

In most cases of our calculations as well as previous 
ones, from the center to the C-0 core edge the oxygen 
abundance has a peak near the outer edge of the oxygen- 
rich core. This configuration is Rayleigh- Taylor unstable 
and will be rehomogenized by convection (Salaris et al. 
1997; Isern et al. 1998). Their results show that the re- 
homogenization proceeds slowly and the effect is relatively 
small during the liquid phase of the WD interior. How- 
ever, because of this effect the final chemical profile of the 
inner SNe la progenitors may differ from the profile at the 
end of second dredge-up depending on the time scale of 
becoming SNe la. 

Figure ^ shows the central C/0 ratio in mass fraction, 
which is not a monotonic function of the metallicity and 
the ZAMS mass. One reason for this complexity is that 
the C/0 ratio is sensitive to the evolutionary change in the 
convective core size during the later stage of He burning. 
Small differences in the core temperature and opacity lead 
to a different evolutionary behavior of the convective core, 
which largely changes the C/0 ratio. 

The central C/0 ratio is sensitive to the ^^C(a,7)^^0 
rate. For example, for M = 5Mq, Z = 0.02, and Y = 
0.27753, this ratio is 0.58 for CF88xl.5 rate, while C/0 = 
0.89, 0.47, and 0.24 for CF88, CF88xl.7, and CF85 rates, 
respectively. This ratio, on the other hand, is not so sen- 
sitive to Y: for example, C/O=0.58 and 0.53 for Y=0.285 
and Y=0.30, respectively. Several groups including Salaris 
et al. (1997) and Bressan et al. (1993) calculated the cen- 
tral C & O mass fractions for solar metallicity. SSMM and 
a series of subsequent papers calculated the fractions for 
various metallicities. Our results of the central C/0 ratio 
is systematically larger than those of SSMM and Salaris et 
al., who adopted the CF85 ^^C(a, 7)^^0 rate and different 
treatments of convection. 
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Figure ^ shows the total mass of ^^C (Af(^^C)) inte- 
grated over the WD mass, Mia — 1-37Mq, just before the 
SN la explosion. Here we assume C/0 '--^ 1 in the outer 
layer of M^o < ^ ^ha. which is formed by He shell 
burning during accretion. This is because our results as 
well as other previous works indicate that a star with C-0 
core mass heavier than ^ IMq tends to yield C/0 ^ 1 
in the product of He shell burning. We note that C/0 
depends on both the (p, T) track and the nuclear reaction 
rate. For example, if we use CF88 for the -^^C(a;, 7)^^0 
rate in the model of Figure 3(b), the C/0 ratio in the 
shell is C/0 ~ 1.46. 

The total carbon mass is important for studying SNe la, 
because it determines the total explosion energy. The total 
carbon mass fraction, M(^^C)/Mia, varies in the range of 
0.36 - 0.5 depending on the ZAMS mass and metallicity, 
which might be related to the observed variation of the 
SNe la brightness. 

In Figures 11 and[l^, we plot the central C/0 ratio and 
Af(^^C) again, but abscissa is changed to Mqq. These 
figures show that by replacing the abscissa from A/ms to 
Mcoi these curves somewhat converge. Its implication is 
discussed in the next section. 

4. SUMMARY AND DISCUSSIONS 

We have calculated the stellar evolution through the car- 
bon ignition or the second dredge-up for M = 3 — QMq 
and Z — 0.001 — 0.03. Using these evolutionary models, we 
have made systematic investigations on the stellar metal- 
licity and mass dependencies of the C-0 core mass Mqq at 
the end of the second dredge- up, the C/0 ratio in the C-0 
core, and the mass range (A/lo, M\jp). The major effect 
of changing metallicity appears in the Mqq — M-c^s relation 
(Figs. 9 and 11). Main conclusions and the implication on 
SNe la are as follows: 

1) The Z dependence of A/up and A/lo a-re important to 
obtain the white dwarf luminosity function, and to model 
the galactic chemical evolution, because these determines 
the SNe la rate. Our results exhibit that Afyp and A/lo 
are both smaller for lower Z . This implies that if the initial 
mass function is independent of metallicity, the number 
density of SN la progenitors are larger in a lower metallic- 
ity environment (Fig. ^); therefore, the luminosity func- 
tion of WDs and SNe la would depend on Z . 

2) Our results indicate that the innermost part of C-0 
WDs is likely to be oxygen-rich (C/0 < 0.6) rather than 
C/0 = 1 as usually assumed in the SNe la modehng. We 
find two tendencies of the central C/0 ratio: (1) For a 
fixed as stellar mass increases, the C/0 ratio at the 
stellar center first increases to the maximum and then de- 
creases. The location of peaks, which depend on metal- 
licity, are seen at lower mass in our results (Mms ~ SM©) 
than in SSMM (M^s ~ TMq). (2) For larger Z, the above 
maximum C/0 ratio is larger. Compared with SSMM, 
our results show a much larger variation in the C/0 ratio 
distribution. 

3) For WDs which can be progenitors of SN la, the central 
ratio varies in the range of C/0 = 0.4 — 0.6 depending 
mainly on the ZAMS mass. The total carbon mass frac- 
tion just before SN la explosion is found to vary in the 
range M{^'^C)/Mi^ = 0.36 - 0.5. These may be partially 
responsible for the observed variations in the brightness 



and spectra of SNe la. 

Most of previous works on the SNe la explosion have 
used progenitor models with C/0=1. Only recently 
Hoflich et al. (1998) have started calculations with dif- 
ferent C/0 for the DD models. The variation of C/0 
ratio would be important, because as the inner regions are 
burned during the deflagration phase, the energy release is 
changed and, consequently, the pre-expansion of the outer 
layers is modified. This changes the total explosion energy, 
and hence the expansion rate of the SNe la, the brightness 
and shape of light curves as well. Also it may affect the 
ptr in the DD model to induce a large variation in the 
explosion energy. 

4) The ranges of the variation in the central C/0 ratio and 
Af(^^C)/Afia obtained above do not significantly depend 
on Z for 0.001 < Z < 0.02. These ratios shown in Figures 
9 and 10 exhibit smaller dispersion with respect to metal- 
licity if these ratios are plotted against Mqq (Figs. 11 
and 12). We conclude that the metallicity effects are rel- 
atively minor for the same Mqo, although there are still 
some significant variations in the central C/0 ratio for 
Mco >0.8Mq (Fig. 11). Moreover, WDs with Z < 0.002 
need not be considered for SNe la progenitors, because 
such low metallicity WDs blow too weak winds to evolve 
into SNe la (Kobayashi et al. 1998). These would be 
good for determining cosmological parameters with the 
use of SNe la, because "evolutionary" effects on the SN 
la properties are minor at different red-shifts. 
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Table 1 

Isotopes included in the nuclear network for He burning. 



Isotope 
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Isotope 
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1 


H 
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He 
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12-13 
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13-15 





14-18 


F 


17-19 


Ne 


18-22 


Na 


21-23 


Mg 


22-27 


Al 


25-29 


Si 


27-30 


Fe 


56 







7 




Fig. 1. — Evolution in the H-R diagram (upper panel) and the central density (pc) vs. central temperature (Tc) diagram (lower panel) of 
our 5 Mq (Z=0.02) models for Y = 0.2775 (solid), compared with SSMM (dotted; Schaller et al. 1992). Dependence on Y is small as seen 
from the case of Y = 0.3 (dashed). A model with a larger mixing length parameter a (=2.25) is shown in a dotted curve with crosses. 
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Fig. 2. — Evolution in the H-R diagrams of the 5Mq models for several metallicities Z (upper panel) and the Z = 0.02 models for several 
stellar masses (lower panel). 
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Flc;. 3. — (a) Abundances in mass fraction in the inner core of the 7 Mq star for Z=0.001 at age 4=48.45 Myr, with the central temperature 
and density log Tc (K)=8.58 and log pc (g cm^'^) =6.82 (left panel). The maximum temperature, log Tmax=8.78, is attained at Mr = 0.54Mq. 
(b) Same as Figure 3(a) but for Z=0.02 at the end of the second dredge-up, t=51.87 Myr, with log Tc=8.29 and log pc =7.04 (right panel). 
Maximum temperature, log rinax=8.69, is attained at Mr = O.72M0. 
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Evolution of 3-9 M„ Stars from Z=0.001 to 0.03 



and metallicity effects on Type la Supernovae 
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M (M^) 

ms ^ 0/ 

Fig. 6. — Metallicity dependence of the ZAMS mass (Mms) vs. C-O core mass (Mco) at the end of the second dredge-up. Filled circles 
indicate that off-center carbon ignition occurs in these model. Hence, the abscissa of these points correspond to the Mup. 
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Fig. 7. — Mup obtained in this work compared with earlier works. The lower ZAMS mass limit for becoming an SN la, Mlqi is also shown. 
Since this value depends on the binary parameters, three curves are shown, which correspond to Mco = 0.75, 0.80 and 0.85 Mq. Crosses 
labeled by 'Sch' are Mup of Schaerer et al. (1993). 
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Evolution of 3-9 M„ Stars from Z=0.001 to 0.03 



and metallicity effects on Type la Supernovae 




Fig. 8. — Metallicity and Mlo dependence of the number of SNe la progenitors normalized to the case with Z = 0.02 and Mlo = O.SMq. 
The Salpeter IMF is assumed. 
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Fig. 9. — Central C/O ratio in mass fraction as a function of metallicity and stellar mass, (left) 

Fig. 10. — Total ^^C mass (M(i^C)) included in the WD of mass Mja = 1.37M0 just before an SN la explosion as a function of M, 
(right) 




